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(57) Abstract: A method for determining component flow 
rates of a multiphase fluid in a conduit (the fluid consisting 
of at least three known components) includes the steps of: 
measuring at each of two different positions along the conduit 
at least four mixture quantities, typically the sound speed, 
the flow velocity of the multiphase fluid, the pressure and 
the temperature; providing a speed of sound in each of the 
components at the measured pressures and temperatures; 
providing a trial value for each of either the component flow 
rates or the phase fractions; using a predetermined model to 
calculate values for the measured mixture quantities based on 
the trial values for each of either the component flow rates 
or the phase fractions; using a predetermined error function 
to determine an error value; and using a predetermined 
optimizing algorithm to determine whether the calculated 
values are acceptable, and, if they are not, to provide a new 
trial value for each of either the component flow rates or the 
phase fractions. In some applications, the error function is the 
sum of the squares of the difference between the measured 
and calculated values at each point. 
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Method and Apparatus for Determining Component 
Flow Rates for a Multiphase Flow" 



The present invention concerns techniques for the measurement of the flow rate 
of each component of a multiphase fluid. More particularly, the present invention 
concerns techniques based on measurements at different locations along a conduit for 
determining the flow rates of each component of a multiphase fluid within the conduit 

In extracting oil and gas from a formation, it is advantageous to monitor the flow 
rates of the different components of the production fluid, usually gas, oil and water. 
Such a three-component mixture presents a problem that is often more difficult to solve 
than that of detennining the flow rates for binary mixtures. For a broad class of well 
mixed binary mixtures flowing in a conduit of known cross sectional area, for which the 
density and the sound speed of each component are known, measuring the sound speed 
of the mixture and the convection velocity of the mixture provides sufficient 
information to determine the flow rates of each of the components. The sound speed of 
the mixture is directly related to the phase fractions of each component in the mixture. 
More specifically, to determine the four unknowns, the phase fraction of each 
component and the flow rate of each component, four equations are used. One equation 
expresses that the sum of the component (phase) fractions is unity. Another is that the 
speed of sound in the mixture depends on the speed of sound and density of the 
individual components and the phase fraction of each constituent. Finally, there are 
two equations relating the flow velocities of the individual components to the overall 
mixture flow rate. 

The approach of using a single sound speed measurement and a mixture flow 
rate measurement used in determining the component velocities in the case, of a binary 
mixture, however, cannot be extended to solve the problem of finding the component 
velocities of a fluid with more than two components. Additional information is required 
beyond what is provided by a measurement of the speed of sound in the mixture and a 
measurement of the flow rate of the mixture. 
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To provide the required additional information, the prior art teaches multiphase 
flow meters that typically rely on several so called orthogonal measurement systems 
(said to be a orthogonal because each measurement system provides information that is 
at least partially independent of the information provided by the other measurement 
systems). Multiphase flow meters according to such an orthogonal system approach 
include meters based on multiple-energy nuclear sources, ultraviolet measurements, 
capacitance measurements, venturi effect measurements, and infrared measurements. 

An alternative approach to determining the component flow velocities of a 
multiphase fluid in a conduit is to determine the additional required information from 
multiple point measurements, i.e. from measurements of the same information at 
different places along the conduit For example, the speed of sound in the fluid, the 
flow rate of the fluid, and the pressure and temperature of the fluid would be made at 
two or more locations along the conduit In addition, a multiphase flow model is used 
to relate the values of the parameters that are measured at one location to those 
measured at another location. These relationships can provide the additional constraints 
required to solve for additional parameters. Typically, the equations are nonlinear in 
several variables. Many methods can be employed to "solve" for the flow parameters of 
interest. 

One class of methods defines an error function based on the -calculated values 
of the parameters compared to the measured values of the parameters. The flow-related 
parameters sought (i.e. the component velocities or component phase fractions) are 
adjusted -iteratively until the error is minimized. In this context, the flow parameters 
that result in minimizing the value of error function are assumed to be correct. 

One company, previously Loke of Norway (now owned by FMOKOS), is 
considered by many to have pioneered the general approach based on multiple point 
measurements for production allocation measurements in oil and gas production 
facilities. One implementation by Loke of the general approach is their software called 
Idun, which uses conventional pressure and temperature measurements along with a 
choke position measurement and knowledge of the fluid property characteristics to 



o 

WO 02/46705 



o 

PCT/GB01/05340 



3 

estimate the component flow rates (or phase fractions). The overall accuracy and 
robustness of the Idun approach is directly influenced by the type and quality of sensors 
available. 

Another approach to determining component flow rates in a multiphase fluid 
based on making measurements at multiple locations along a conduit carrying the 
multiphase fluid is that based on a gradio-venturi system, which includes a venturi 
meter and employs a remote pressure sensor located several hundred feet above the 
venturi. The pressure difference between the pressure at the venturi and that at the 
remote transducer can be related to the flow rate and composition through a multiphase 
flow model, and can be used in conjunction with the pressure difference due to the flow 
thorough the venturi to estimate the component flow rates. Such an approach has 
several drawbacks. It requires a venturi, which is intrusive, and it has an accuracy 
limited in two-phase flow to ±10% of the total flow. Moreover, the accuracy degrades 
substantially in the presence of any significant entrained gas. 

What is needed is a system of measuring the component velocities and phase 
fractions of a fluid that includes at least three components, and that is nonintrusive, 
sufficiently accurate, and that does not provide spurious solutions because of 
insufficient information. 

Accordingly, the present invention provides an apparatus and corresponding 
method for determining component flow rates of a multiphase fluid in a conduit, the 
fluid consisting of at least three known components, the method including the steps of: 
measuring at each of two different positions along the conduit at least four mixture 
quantities; providing a speed of sound in each of the components at the measured 
pressures and temperatures; providing a trial value for each of either the component 
flow rates or the phase fractions; using a predetermined model to calculate values for 
the measured mixture quantities based on the trial values for each of either the 
component flow rates or the phase fractions; using a predetermined enror function to 
determine an error value; and using a predetermined optimizing algorithm to determine 
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whether the calculated values are acceptable, and, if they are not, to provide a new trial 
value for each of either the component flow rates or the phase fractions. 

In a further aspect of the invention, the error function is the sum of the squares 
of the difference between the measured and calculated values at each point. 

In another, further aspect of the invention, the four mixture quantities are the 
sound speed, the flow velocity of the multiphase fluid, the pressure and the temperature. 

The above and other objects, features and advantages of the invention will 
become apparent from a consideration of the subsequent detailed description of a 
preferred embodiment of the invention given, by way of example, with reference to the 
accompanying drawings, in which: 

Figure 1 is a schematic/ block diagram of a system for measuring component 
flow velocities for a multiphase fluid in a conduit, according to the present invention; 

Figure 2 is a set of plots showing the results of a simulation of use of the present 
invention in case of no free gas being present at an upstream position on a conduit; 

Figure 3 is a set of plots showing the results of a simulation of use of the present 
invention in case of 20% (a phase fraction) free gas being present at an upstream 
- position on a conduit; and - . 

Figure 4 is a set of plots showing the results of a simulation of use of the present 
invention in case of 40% free gas being present at an upstream position on a conduit 

The present invention for determining component flow rates (or, equivalently, 
phase fractions) for a multiphase fluid in a conduit is based on making measurements at 
different points along the conduit, as indicated in Figure 1. The term flow rate is used 
here in reference to a point along a conduit to mean mass per unit time flowing past the 
point along the conduit The flow rate at a point is related to the average flow velocity 



o 

WO 02/46705 



o 

PCT/GB01/05340 



5 

at the point by the cross sectional area of the conduit at the point. The term phase 
fraction of a component of a multiphase fluid in a conduit (at a point along the conduit) 
is used here to indicate the fraction of the total mass of a sample of the fluid (at the 
point along the conduit) that is due to the presence of the component. Knowing the 
flow rate of each component is equivalent to knowing the phase fractions and the 
mixture flow rate. The flow rate for a component, at a point along a conduit, is simply 
the product of the component flow velocity, the phase fraction for the component, and 
the total mass per unit length at the point along the conduit. 

In the particular embodiment of a system according to the invention indicated in 
Figure 1, four mixture quantities, namely the sound speed, the (mixture) flow velocity, 
the pressure and the temperature, are measured at two different positions along the 
conduit 1 1, at position A and at position B, using a measurement module 12 at A and a 
measurement module 13 at B, respectively. In the preferred embodiment, it is assumed 
that the flow can be modeled as time stationary so that measurements at the two 
positions at the same time are considered to be related. In other embodiments, however, 
the assumption that the flow is time stationary can be relaxed. In such embodiments, 
measurements at position B (downstream of position A) would be assumed to be related 
to measurements at position A only after a time interval depending on the flow velocity 
measured at position A, and possibly also at position B. 

Still referring to Figure 1, according to the invention values of the measured 
quantities are provided to a multiphase flow modeler 15, which also receives _as input 
the speed of sound in each of the individual components. Such component speed of 
sound values are determined based on pressure and temperature measurements at 
positions A and B and on an assumed PVT model. The multiphase flow modeler 15 
then relates the measurements at position B to the measurements at position A, 
according to the methodology described below, which involves using successive trial 
values of the component flow rates (or the phase fractions) and an error function that 
compares the measured mixture quantities with calculated mixture quantities based on 
the trial values, and iterating until the value of the error function is acceptably small. In 
the preferred embodiment of the invention, the system includes an error function 
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evaluator 16 that computes, for each new trial value of the component quantities 
(component flow rates or phase fractions), each new corresponding value of the error 
function. An optimizer 17 then determines whether the calculated values are close 
enough to the measured values based on a predetermined tolerance, and if not, 
determines new trial values for the component quantities* The component rates (or the 
phase fractions) that minimize the error are ultimately provided by the optimizer 17 as 
the output of the system. 

For an N-component mixture in a conduit, the component densities p t , phase 
fractions $ t , and component sound speeds c k are related to the mixture sound speed 
and mixture density by the equations, 

N 



Pmix = 2 tiPi < CD 



1=1 

_J_- V * 

2 - 2Li „Ji . and (2) 

Pmvc C mix 1=1 PA 

2>< =1 . <3> 

For fluids in a flexible conduit, the propagation velocity of sound, or any other 
compressional wave, is influenced by the structural properties of the conduit. For a 
fluid in a conduit that is a vacuum-backed tube of radius R , thickness t and Young's 
modulus E 9 the measured mixture sound speed c meas is related to the infinite domain 
sound speed c mix and to the structural properties of the tube via the relation, 

1 1 2R 

~~2 =2 . (4) 

However, as mentioned above, for a mixture of three or more components (called here a 
multi-component mixture, as opposed to a binary component mixture), the speed of 
sound no longer uniquely identifies the phase fractions of the different components. 
Additional information is required. 
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The above relations combine to provide a relation between measured sound 

2R 

speed c meas , input fluid and structural properties represented by & — , and unknown 

phase fractions p t (where the index refers to one or another of the components), a 
relation that can be expressed in the form, 
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where the index a refers to either position A or position B, and the index i refers to the 
different components of the fluid. 

The invention will now be described in an embodiment that is particularly 
advantageous for determining the component flow rates (or, equivalently, the phase 
fractions) for a production fluid being pumped from a formation containing oil, gas, and 
water. It should however be understood that the invention has more general 
applicability; as will be clear, nothing about the invention restricts it to fluids of oil, gas, 
and water. 

In the particular case of a three-phase (three-component) fluid consisting of gas, 
oil, and water, equation (5) becomes, 

1 



2 rmix,a 
mix y a 



T + T + ~ + a a 

Pm C ga Poa C oa Pwa C y 



ga ga r oa oa r^wa^wa 

with <P za + 0oa + 0wa = ^ > ^ where the index a refers to either position A or to 
position B along a conduit bearing the production fluid. 

The invention imposes, as one condition, mass flow continuity of the individual 
phases, and uses a model for the evolution of gas from the oil. As the production fluid 
is pumped from position A up to position B, closer to the surface, gas escapes out of the 
oil/gas solution. The present invention uses a model that calculates the amount of gas at 
position B based on, among other factors, the amount of gas and oil at position A. 
Conceptually, the model tracks the multiphase flow as a four-component mixture, a 
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mixture consisting of water, oil, gas present at position A, and gas evolved from the oil 
between position A and position B . 

The invention will now be described in the particular case of well-mixed flow of 
a fluid consisting of oil, gas, and water. The extension to the general case, i.e. to the 
case of flow that is not well-mixed, will be apparent to one skilled in the art. 



In the case of well-mixed flow, there is no slip among the individual phases, i.e. 
the mixture is homogeneous. In such a mixture, the condition of mass flow continuity 
from position A to position B becomes, 

Pmix i A^A U A = Pmix 9 B^B U B , ( 7 ) 

where A a is the cross section of the conduit at the position indicated by the subscript a , 

and u a is the corresponding flow rate. In terms of the individual phases, assuming that 

the cross sections at position A and position B are equal, the condition of mass flow 
continuity expressed by equation (7) becomes, 

U A^oaPoA + 0gAPgA + ^ W APwAh^oBPoB + ^ g 3PgB +0wbPwb) • (8) 

Conservation of the individual phases results in three equations, one equation for 
each of the three components. For the oil and water components, the invention in the 
preferred embodiment assumes that the density of the water and oil components is the 
same at position A and position B, yielding, 

U A0oAPoA=U B 0 oB P oB , (9) 

and 

u a0«>aPwA = u b0wbPwB . ( 10 ) 

These simplifying assumptions can be relaxed without altering the basis of the approach 
of the present invention. 



Hie third equation derives from the conservation of the gas phase. The present 
invention models two mechanisms to account for the change in volumetric flow rate of 
the gas. First, it is assumed that the gas present at position A under goes an expansion 
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characterized by its compressibility factor. Secondly, it is assumed that amount of gas 
that evolves from the fluid as it moves from position A to position B is proportional to 
the amount of oil present at position A. Expressing these relations in terms of 
volumetric flow rates at standard temperature and pressure (indicated by the subscript 
R)_ yields the following relations, 

**A 2 A B r R J B 

where f} is the so-called gas evolution constant, which is defined as the amount of gas 
(in standard cubic feet) evolved from a barrel of oil as the pressure is lowered one 
pound per square inch (psi), and Z is the gas compressibility factor, defined as, 

V R P R T 

The gas evolution constant fi can be determined from a standard differential 
vaporization test in which the amount of dissolved gas (in standard cubic feet) per barrel 
of residual oil is reported as a function of pressure (see e.g. PVT and Phase Behavior of 
Petroleum Reservoir Fluids, by Ah Danesh, Elsevier, 1998). The gas compressibility 
factor Z is also reported in typical PVT field analyses of black oils. 



The first term in equation (11) represents the amount of gas present at position 
A, and the second term represents the amount of gas evolved from the oil as the fluid 
moves from position A to position B. 

The invention therefore constructs a system of seven equations. The sound 
speed mixing laws, given by equation (6), provide two equations. The phase 
conservation conditions, given by equations (9-11), provide three equations. These 
seven independent equations are in terms of six unknowns (the component flow rates or 
the phase fractions at both positions A and B) and so provide more than enough 
information to uniquely determine the six unknowns. Since the system of equations is 
over determined, instead of attempting to solve the system of equations (which may 
include inconsistent conditions), the invention seeks a solution that minimizes an error 
function. In the preferred embodiment, the invention uses as an error function a simple 
sum of the squares of the differences between the calculated values of a quantity 
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measured at the two positions A and B and the measured values, i.e. using F t to 
represent a quantity that is measured at position A and position B, where the index i 
ranges over N such quantities, the invention uses as the error function, 

E-j^F^-F^l (13) 

The application of this general approach can take many forms, i.e. many different 
algorithms or procedures may be used to minimize the error function E; how, 
specifically, the error function is minimized is not the subject of the present invention. 

The accuracy of the general approach used in the present invention, that of 
minimizing an error function as explained above, depends on many factors, including: 
the fidelity of the parametric model; the quality, quantity and types of measurements; 
the accuracy of the calibration of the model input parameters; the accuracy of the well 
parameters, including for example the radius R , thickness * and Young's modulus E 
of the conduit; and the accuracy of the fluid parameters, including for example the 
component speeds of sound 

The inventors conducted several Monte Carlo simulations to evaluate the 
multiphase flow measurement approach developed above. The well parameters were 
selected to be representative of a producing black oiL The parameters used in the 
simulation are listed below. 

Table 1. Parameters of components. 



Fluid 


Density 
(kg/m*3) 


sound speed 
(ft/sec) 


oil 


700 


3500 


water 


1000 


5000 


gas 


l^P/Pata 


1200 
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Table 2. Parameters of the mixture and conduit 



Component-independent parameters 


Pressure at A 


3000 psi 


Pressure at B 


2800 psi 


Temperature at A and B 


100 C 


P, gas evolution constant 


0.316 scffstb/psi 


Inner diameter of conduit 


3.0 inches 


Wall thickness at A 


0.50 inches 


Wall thickness at B 


0.25 inches 



Figures 2-4 present the results of the Monte Carlo simulations, in which the 
input parameters and measured parameters are varied in a random, but bounded, 
manner. The three figures illustrate the ability of the invention, even using the 
simplifying assumptions noted above, to measure the multiphase flow parameters as a 
function of oil and water cut for three levels of free gas present at position A (namely 
0% for Figure 2, 20% for Figure 3, and 40% for Figure 4) in the presence of the 
uncertainty for each of the parameters indicated in Table 3 below. The pressure was not 
directly varied in the simulations; the effect of any uncertainty in pressure is accounted 
for in the variations of the other model parameters. 

Table 3. Uncertainty in model parameters. 



Parameter 


Uncertainty 


P, gas evolution constant 


2% 


Oil sound speed 


2% 


Gas sound speed 


1% 


Water sound speed 


1% 


Oil density 


1% 


Gas density 


2% 


Water density 


0.5% 


Measured sound speed at A and B 


0.25% 


Measured flow velocity at A and B 


0.25% 
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For each of the simulations indicated in Figures 2-4, a single wall sub (measurement 
station) is assumed, with a thickness at position A of 0.5", and a thickness at position B 
of 0.25", as indicated in Table 2. 

The simulations indicate that the invention provides reasonable estimates of oil, 
gas, and water phase fractions, within ±10% of reference values for a reasonable level 
of uncertainty in both measured and input parameters for free gas volume fractions 
ranging from 0% at position A to 40% at position A 

The solution appears to be robust to uncertainty; gradually increasing the 
uncertainty in the input and measured parameters appears to translate into only a 
gradual degrading in flow measurement accuracy. 

It is to be understood that the above-described arrangements are only illustrative 
of the application of the principles of the present invention. Numerous other 
modifications and alternative arrangements may be devised by those skilled in the art 
without departing from the spirit and scope of the present invention, and the appended 
claims are intended to cover such modifications and arrangements. 
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CLAIMS: 

1. A method for determining component flow rates of a multiphase fluid in a 
conduit, the fluid consisting of at least three known components, the method comprising 
the steps of: 

a) measuring at each of two different positions along the conduit at least four 
mixture quantities; 

b) providing a speed of sound in each of the components at the measured 
pressures and temperatures; 

c) providing a trial value for each of either the component flow rates or the phase 
fractions; 

d) using a predetermined model to calculate values for the measured mixture 
quantities based on the trial values for each of either the component flow rates or the 
phase fractions; 

e) using a predetermined error function to determine an error value; and 

f) using a predetermined optimizing algorithm to determine whether the 
calculated values are acceptable, and, if they are not, to provide a new trial value for 
each of either the component flow rates or the phase fractions. 

2. A method as in claim 1, wherein the error function is the stun of the squares of 
the difference between the measured and calculated values at each point. 

3. A method as in claim 1 or 2, wherein the four mixture quantities are the sound 
speed, the flow velocity of the multiphase fluid, the pressure and the temperature. 

4. A method for determining component flow rates of a multiphase fluid in a 
conduit, the fluid consisting of at least three known components, the method comprising 
the steps of: 



g) measuring at each of two different positions along the conduit at least four 
mixture quantities; 
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h) providing a speed of sound in each of the components at the measured 
pressures and temperatures; 

i) providing a trial value for each of either the component flow rates or the phase 
fractions; 

j) using a predetermined model to calculate values for the measured mixture 
quantities based on the trial values for each of either the component flow rates or the 
phase fractions; 

k) using a predetermined error function to determine an error value; and 
1) using a predetermined optimizing algorithm to determine whether the 

calculated values are acceptable, and, if they are not, to provide a new trial value for 

each of either the component flow rates or the phase fractions. 

5, A method as in claim 4, wherein the error function is the sum of the squares of 
the difference between the measured and calculated values at each point 

6. A method as in claim 4 or 5, wherein the four mixture quantities are the sound 
speed, the flow velocity of the multiphase fluid, the pressure and the temperature. 
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